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1. Introduction 

The title of this talk refers to wee par- 
tons in large nuclei. The first part of this 
section will be a brief discussion of the the- 
ory and phenomenology of wee parton dis- 
tributions in QCD. The latter part of this 
section (and subsequent sections) will be 
about wee partons distributions in large 
nuclei and their relevance to experiments 
at RHIC and LHC. 

1.1. A brief introduction to the low 
x problem in QCD 

One of the more interesting problems in 
QCD today is to understand the space- 
time (or equivalently, momentum) distri- 
butions of "wee partons" in hadrons and 
nuclei. The phrase "wee parton" was 
coined by Richard Feynman to de- 
scribe those partons which carry a very 
small fraction of the longitudinal momen- 
tum of the parent hadron or nucleus. This 
momentum fraction is the Feynman x 
variable and for wee partons, x<l. To 
lowest order in the coupling constant, the 
Feynman x variable is equal to the well 



known Bjorken x of deeply inelastic scat- 
tering experiments. 

Why is this problem of wee par- 
ton distributions interesting? To un- 
derstand this let us first briefly review 
the usual Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi branching process in per- 
turbative QCD ||. This branching pro- 
cess describes a virtual, ladder-like cas- 
cade, where partons with higher values of 
transverse momentum k t and x, "split" 
into softer partons with lower values of k t 
and x and so on along the ladder. Each 
rung in the ladder provides logarithmic 
contributions to the cross section of the 
sort 



as 



d 2 k t 



dx 

x 



(1) 



In the Double Leading Logarithmic Ap- 
proximation (DLLA), both x and k t are 
ordered along the ladder; this is equivalent 
to summing leading logs in both variables. 
In the DLLA, the gluon distribution func- 
tion, for x <C 1, grows as |J 



dN I ^2 „ 

ax 



exp 2Vlog(l/x)loglog(Q 2 /Q§) , (2) 



*Plcnary talk at Quark Matter '95, 
CA. 



Monterey, 



where Ql and Q 2 are the respectively the 
lowest and highest values of k 2 along the 
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rungs of the ladder. The latter value 
is usually identified with the momentum 
transfer squared in deeply inelastic scat- 
tering experiments while the former typi- 
cally sets the factorization scale. 

However, if x is very small and Q 2 is 
not large enough, the leading logarithmic 
terms in x will dominate and it may be 
appropriate to sum those alone. In this 
approximation, known as the BFKL ap- 
proximation [Q (after the initials of the 
authors), one has ordering along the lad- 
der rungs only in x. In the Regge-Mueller 
language, this corresponds to a t-channel 
exchange of a "perturbative pomeron" in 
the virtual photon-l- target scattering pro- 
cess. The gluon distribution then has the 
form 

dN 2 2 [Q* 1 
x— (x,Q ;Qo)~\Jq2^, (3) 

where A = 4a s log(2). This behaviour of 
the gluon distributions is often called the 
Lipatov enhancement. An interesting fea- 
ture about the BFKL approach is that, 
unlike the DLLA, it is not a twist expan- 
sion (as in the Operator Product Expan- 
sion). 

Regardless of whether the DLLA or the 
BFKL approximations are more appro- 
priate for the x values of interest, they 
have one outstanding feature in common: 
for i < 1, the density of wee partons 
(dN/dx) grows rapidly. In both cases, 
it grows faster than 1/x. At sufficiently 
small values of x, one enters a regime 
where partons from neighbouring ladders 
overlap spatially. This regime of "high 
density" QCD is of interest because it 



describes the many body behaviour of the 
quanta of the fundamental theory. 

One consequence of the overcrowding of 
partons is that two soft partons may re- 
combine to form a harder parton. Another 
consequence of the overcrowding is the 
screening of parton-parton interactions by 
the cloud of surrounding wee partons [|] . 
Both of these processes may inhibit the 
growth of parton distributions and cause 
them to saturate at some critical value 
of x. Incidentally, the saturation of par- 
ton distributions clS X — > is required by 
unitarity-the so-called Froissart bound. 
Because these many body phenomena are 
very complex, it is especially gratifying to 
realize that they may (to a degree) be un- 
derstood using well known, weak coupling 
techniques . 

Why are weak coupling techniques ap- 
plicable in this problem? The only scale 
in the problem at low x is p, the den- 
sity of partons. The coupling constant as 
will run as a function of this scale and if 
Aqcd <S Pi then cxs(p) ^ 1- This qual- 
itative argument can be proven rigorously 
and is the basis of all that follows. 

Thus far, we have focused on some 
of the theoretical reasons why one must 
better understand wee parton distribu- 
tions. These theoretical impulses have 
been around from the pre-history of QCD. 
What is re-generating interest in a previ- 
ously moribund field is a new generation of 
experiments, which have begun to probe 
the high density phase of QCD. Fore- 
most among these is HERA, the electron- 
proton collider experiment at DESY ||. 
At HERA, x's as low as x ~ 10~ 5 and 
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Q 2 as high as Q 2 ~ 10 5 GeV 2 are at- 
tained. One of the observables measured 
at HERA is F% p (x, Q 2 ), the electromag- 
netic form factor of the proton. At the 
above mentioned values of x, F% continues 
to rise with decreasing x, for a wide range 
of Q 2 values. 

Since F 2 can be related directly to the 
sea quark distributions and indirectly to 
the gluon distributions, these experimen- 
tal results seem to confirm the behaviour 
predicted by the theoretical models dis- 
cussed above. The data are apparently 
not good enough to distinguish between 
the DLLA and BFKL scenarios-this issue 
is controversial and will hopefully be set- 
tled conclusively in the near future. 

Another class of experiments which 
probe the region of low x are the deeply in- 
elastic scattering experiments at Fermilab 
(E665) and CERN (NMC). These experi- 
ments are fixed target experiments where 
a high energy lepton scatters off a nuclear 
target. It is observed || that the ratio 
F^/AF? shows a significant depletion at 
x < 10~ 2 - a phenomenon known as "shad- 
owing". The shadowing of nuclear distri- 
butions at low x shows an unambiguous 
A dependence of A ' 9 . Can shadowing be 
understood as a consequence of the faster 
(i.e., at higher x) saturation of nuclear dis- 
tributions as opposed to hadron distribu- 
tions? Naively, the answer is yes, but a 
conclusive answer will have to be a quan- 
titative one. 

Finally, there are the collider experi- 
ments at RHIC and LHC which will also 
probe fairly low x's (~ Pt/y/s). The wee 
parton distributions of the nuclei will pro- 



vide the initial conditions for the evolution 
of the system into a quark gluon plasma. 
Further, most of the observables in a nu- 
clear collision will depend sensitively on 
the wee parton distributions. In the fol- 
lowing sub-section we will discuss briefly 
the problem of wee parton distributions 
in large nuclei. Subsequent sections will 
address the problem in greater detail. 

1.2. Wee partons in large nuclei: an 
outline 



In recent papers ||10| , pT| , Larry McLer- 
ran and I have argued that large nuclei are 
an excellent theoretical tool to study the 
low x problem. We show that for large 
nuclei, in the range of x <C A^ 1 ^ 3 , the 
large density of partons introduces a scale 
/i 2 « A 1 / 3 fm~ 2 , which can be identified 
as the average valence quark color charge 
squared per unit area. When [x 2 is large 
(A is large) the coupling constant «s(/U 2 ) 
is small and a weak coupling expansion 
can be employed to compute quark and 
gluon distributions. 

On the phenomenological side, we ex- 
pect that this first principles calculation 
of quark and gluon distributions in QCD 
will help clarify our understanding of the 
initial stages of nuclear collisions at RHIC 
and LHC. There exist state of the art 
partonic cascade models which provide an 
excellent initial orientation to this prob- 
lem ||12| . However, since cascade models 
do not accurately simulate the wee par- 
ton sea, the enhancement of the momen- 
tum scale due to the large parton den- 
sity will, in comparison, likely lead to i) 
higher temperatures, ii) enhanced minijet 
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production and iii) an enhanced contribu- 
tion of "intrinsic" strangeness and charm, 
than those predicted in these models. 

We will now discuss how the paper is 
organized. We formulate the problem of 
low x parton distributions as a many body 
problem in the presence of an external 
source. The source here is a sheet of static, 
valence quarks localized in the longitudi- 
nal direction but infinite in extent and 
uniformly varying in the transverse direc- 
tion. Next, we write down a partition 
function for this system which includes 
a stochastic averaging over the external 
sources of color charge. This procedure 
introduces an additional dimensionful pa- 
rameter fi 2 in the theory. It is defined as 
the average valence quark colour charge 
squared per unit area. 

The approach we adopt to solve the 
above mentioned problem is as follows: we 
first obtain the classical background field 
in the presence of the source. It turns out 
that there is a very simple solution to the 
classical field equations in which the only 
non-zero fields are the transverse vector 
fields which are zero ahead of the source 
and two dimensional pure gauge fields be- 
hind the source. The gluon distribution 
functions in the background field can thus 
be expressed as correlation functions of a 
two dimensional, Euclidean field theory. 

For parton 
transverse momenta a 2 fi 2 << k 2 , which 
we identify weak coupling regime, the 
gluon distributions obey the well known 
Weizsacker-Williams distribution scaled 
by jJL 2 ~ A 1 / 3 . In the strong coupling 
regime k 2 << a 2 [i 2 , the distribution func- 



tions still vary as 1/x but the transverse 
momentum dependence is changed. This 
likely reflects the sharp decay of correla- 
tion functions in transverse space. We 
have not succeeded in obtaining an ana- 
lytical expression for the transverse mo- 
mentum dependence-but have succeeded 
instead in writing an algorithm to solve 
the equations numerically using a Monte 
Carlo procedure. 

We next outline the procedure to com- 
pute small fluctuation Green's functions 
for scalars, vectors and fermions in this 
classical background field. The detailed 
computations, performed in collaboration 
with A. Ayala, J. Jalilian- Marian and L. 
McLerran are reported elsewhere [Hi|^H> . 
These Green's functions are useful when 
one computes the higher order contribu- 
tions in as to the quark and gluon dis- 
tribution functions. The fermion Green's 
functions are used to determine the sea 
quark distributions. We find that the 
enhanced momentum scale arising from 
the higher density of partons leads to an 
enhancement of the strange and charm 
quark contribution to the nuclear wave- 
function. 

We then proceed to discuss briefly 
higher order corrections to our results - 
specifically, the one loop corrections to the 
classical background field and the origin of 
the Lipatov enhancement in our approach. 
We also speculate how our solutions may 
be iterated to all orders in as- Finally, we 
discuss recent progress in extending our 
approach to nuclear collisions and end our 
presentation with a brief summary. 



5 



2. Formulation of the problem 

In this section we will formulate the 
problem of calculating the distributions of 
partons in the nuclear wavefunction 
well defined many body problem. To do 
this we will work in the Infinite Momen- 
tum frame and use the technique of Light 
Cone quantization. We will also be work- 
ing in Light Cone gauge A + = 0. There 
are several advantages to these choices. 
One of these is the fact that it is only 
in Light Cone gauge that partons have 
a manifest physical interpretation as the 
quanta of the theory. Further, if we work 
on the Light Cone, it is possible to con- 
struct a simple and intuitive Fock space 
basis [16| on which our formalism relies 
heavily. One consequence is that the elec- 
tromagnetic form factor of the hadron 
which is measured in deeply inelastic scat- 
tering experiments, is simply related to 
parton distributions by the relation [17 



the scope of this talk. For an alternative 
approach to the problem, the reader may 
consult the recent paper by Makhlin 



2.1. A static source on the Light 
Cone 

We begin by making several physically 
plausible assumptions. First, if the den- 
sity of partons per unit area at low x 
is large, and weak coupling techniques 
are applicable, the recoil experienced by 
the valence quarks due to bremsstrahlung- 
ing of low x partons is small-the valence 
quarks obey straight line trajectories. We 
can thus replace the valence quarks in the 
infinite momentum frame nuclear wave- 
function by static, external sources of 
charge. Next, if we require that the nu- 
cleus be Lorentz contracted to a size which 
is much smaller than the wavelength of the 
parton in a frame co-moving with its lon- 
gitudinal momentum, then 



2Rm/P « 1/xP x « 1/Rm, (5) 

2 2 + k + + 

F 2 {x,Q ) = (J d k t dk + x5(x - — ) ^ a A a A ) w i|4) e ^> j s t h e nuc i ear ra dius, P the mo- 

w_ mentum of the nucleus and m the nu- 
cleon mass. From the above inequality, 
we deduce that when x << A^ 1 ^ 3 , the 
partons see a sheet of static color charges 
of large transverse extent but localized in 
the longitudinal direction. Since the av- 
erage transverse momentum scale k± » 
l/R, the valence quark distribution can be 
taken to be uniform and infinite in extent 
in the transverse direction. 

It is appropriate to use Light Cone co- 
ordinates to investigate the dynamics of 
the infinite momentum frame wavefunc- 
tion of the nucleus. The conventions used 



In the above, P + is the momentum of the 
nucleus, k + = (k + k z )/y2 and k t are 
the parton "longitudinal momentum" and 
transverse momenta respectively, Q 2 is the 
momentum transfer squared from the pro- 
jectile and a^a is the number density of 
partons in momentum space. Hence, in 
principle one only needs up to integrate 
the calculated distributions to the scale 
Q 2 to make comparison with experiment. 
One difficulty in the Light Cone formal- 
ism is the problem of choosing appropriate 
boundary conditions. This issue is beyond 
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here are discussed in Ref. lOfl. In the 



Light Cone formalism, the "external" cur- 
rent experienced by the low x partons due 
to the static, valence quark charges takes 
the form 



the ground state of the low x partons in 
the presence of the valence quark external 



source is 



JZ = ^ + Pa(x + ,X ± )5(x-). 



(6) 



In the above, x ± = (t ± x)J y/2, where x + 
is the Light Cone time. 

The Light Cone Hamiltonian P ( 
which generates translations in x + ), in 
the presence of an external source, can be 
written as 

P- = \F t 2 + ^(p F + D t .E t )^(p F + D t 

Here pp is the charge density due to 
the external source plus the dynamical 
quarks, A t and tp are the dynamical vector 
and spinor fields, E t = d-A t is the trans- 
verse electric field and F 2 is the trans- 
verse field strength tensor. In writing the 
above expression, we have made use of the 
constraint equations on the Light Cone to 
eliminate the non-dynamical fields. 

2.2. Ground state expectation val- 
ues in the presence of external 
sources 

In QED, the infinite momentum frame 
wavefunction of the system with the ex- 
ternal source in Eq. (|^) is a coherent 
state IfLOl. We have not succeeded in do- 



ing the same in QCD and have concen- 
trated instead on computing ground state 
expectation values in the presence of ex- 
ternal sources. The partition function for 



Z=<0\e iTP |0>= lim V < N\e iTP ~ \N > Q . (8) 



N 



The sum above also includes a sum over 
the color labels of the source of color 
charge (denoted by Q) generated by the 
valence quarks. 

In principle, for quantized sources of 
color charge, evaluating the trace over dif- 
ferent values of the color charge is diffi- 
cult. However, for a large nucleus, the 
problem can be simplified considerably. 
Visualize the transverse space as a grid 
of boxes of size d 2 x t » l/p va i e „ C e ~ 
j$7) or equivalently, of parton transverse 
momenta q 2 « A 1 ^ MeV 2 . In this 
kinematic range, the number of valence 
quarks, and therefore the charge Q in each 
box, will be much greater than one. A 
large number of charges corresponds to a 
higher dimensional representation of the 
color algebra-the sum of the color charges 
of the valence quarks in a grid can be 
treated classically. 

Further, if the total charge in the grid 
is much less than the maximum possible 
charge the central limit theorem tells us 
that the density of states corresponding 
to charge Q is a Gaussian e~® / 2m . Sum- 
ming over the color labels of the states is 
therefore equivalent to introducing in the 
path integral the measure 

[dp}exp(-^ [ d 2 x t p 2 (x t )) . (9) 



2p 2 „ 
This ensures that 

< p a {x t )p\y t ) >= p 2 5 ab 5®(x t - y t ) ; < p a (x t ) >= 
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Above, the variance /i 2 is the average va- 
lence quark color charge squared per unit 
area. It can be written as 

3A 4 
iiR 2 3 

where < Q 2 >= 4g 2 /3 is the average 
charge squared of a quark. 

We can now write the path integral rep- 
resentation of the partition function Z in 
the Light Cone gauge A__ = as 

Z = J [dA t dA + ][d^dip}[dp\ 

exp [iS + ig / d 4 xA + (x)5(x~)p(x) — - 



2. 



Hence, the result of our manipulations 
is to introduce a dimensionful parameter 
p? « 1.1 A 1 / 3 fm -2 in the theory. If we im- 
pose the current conservation -D M J M = 
and integrate over the external sources p, 



we obtain an expression JT0[ containing 
modified propagators and vertices. 

We shall however attempt to solve the 
above many body problem in the follow- 
ing fashion: we first solve for the classical 
background field, compute the small fluc- 
tuation determinant in the background 
field to obtain Green's functions and prop- 
agators, use these to compute the first ra- 
diative corrections to the background field 
and systematically iterate the process to 
all orders. 

3. The Classical Background Field 

We first discuss the problem of comput- 
ing the solution of the classical equations 
of motion for the gluon field in the pres- 
ence of a source which is a delta function 



along the Light Cone. The second part 
of this discussion is about computing cor- 
relation funtions in this background field. 
Notice from Eq. ([12]) that there is still an 
up^gral over p with a Gaussian weight to 
average over at the end. At each order in 
as, this procedure is equivalent to sum- 
ming the theory to all orders in the effec- 
tive coupling dsn/kt- 

3.1. Solutions to the Classical Equa- 
tions 



The equations of motion are 



(13) 



where J, the classical Light Cone source, 
may be represented as 



JZ = S + »Q a (x + ,x t )5(x- 



(14) 



We will work in Light Cone gauge where 
A_ = -A + = 0. 

There exists a solution of the equations 
of motion for this problem, where the lon- 
gitudinal component A + , which is not zero 
by a gauge, vanishes by the equations of 
motion 



A+ = -A' = . 



(15) 



The only non-zero components of the field 
strength are the transverse components 
which we require to be of the form 



Ai(x) = 9(x )ai(x t ) 



(16) 



This functional form is what we expect for 
a classical field generated by a source trav- 
eling close to the speed of light with x = t. 
For x > t, the source has not yet arrived, 
and for x < t, the source should produce 
a field. 
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If we further require that F %3 = 
(where i and j are transverse components) , 
we see that we have a solution of the equa- 
tions of motion as long as 



V • a = gp{x t ) . 



(17) 



Here p is the surface charge density asso- 
ciated with the current J. There is no de- 
pendence on x~ because we have factored 
out the delta function. The dependence 
on x + goes away because of the extended 
current conservation law D^J^ = in the 
background field. 

The condition that i™ = is precisely 
the condition that the field a is a gauge 
transform of the vacuum field configura- 
tion for a two dimensional gauge theory. 
The field configuration which is a gauge 
transform of the vacuum field configura- 
tion for a two dimensional field theory 
may be written as 

t-oh = -—UViU ] . (18) 
W 

The dependence of the U fields on the sur- 
face charge density p is then given by the 
relation 



V • UVU ] = -ig 2 p(x t ) 



(19) 



We have not been able to construct ana- 
lytical solutions for the above equation for 
arbitrary dependence of the surface charge 
density on xt- However, recently we have 
found a numerical algorithm which will 
solve the above equation on a 2-D lattice 
using a Monte Carlo procedure [[Hfl. 

3.2. Computing Correlation Func- 
tions 

To compute correlation functions as- 
sociated with our classical solutions, we 



must integrate over all color orientations 
of the external sheet of charge. This is 
equivalent to computing the expectation 
value of (ai(x t )aj(y t )) with the measure 
(expressed in terms of the compact fields 
U) 



[dU] exp ( j-^ / d 2 x t Tr 



v-(c/ T w f ) 



We see that the measure for this theory 
is that for a two dimensional Euclidean 
field theory where the expansion param- 
eter is a 2 s p 2 /k 2 . Since the theory is two 
dimensional, it should be straightforward 
to compute correlation functions to all or- 
ders in a 2 s p 2 /k 2 using lattice Monte Carlo 
methods. Further, since the theory is ul- 
traviolet finite, there should be no prob- 
lems extrapolating to the continuum limit. 

Since we can solve Eq. ([T9|) numeri- 



cally, the correlation functions can still be 
expressed as an integral over the Gaus- 
sian measure in p. We then do not need 
to worry about the nasty Fadeev-Popov 
term in the above measure. 

The correlation function for the com- 
putation of the transverse momentum de- 
pendence of the structure function may be 
simply evaluated to be 



D{k t ) = —H{k t ) 
as 



(21) 



where D is the propagator for the two di- 
mensional theory 



(2tt) V(fc t - q t )5 tJ D(k 



< aiOj > (22) 



The relation between distribution func- 
tions and propagators is straightforward 
and is discussed explicitly in Ref. (13). For 
instance, for a scalar field <p a ( 



det(V-i 
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the fundamental representation, the dis- 
tribution function is given by the relation 



dN 



2 1 



. (2k- 



(2vr) : 



J2D aa (x + ,k, x + ,k).(23) 



Correlation functions of x t should die ex- 
ponentially at large distances, or alterna- 
tively the Fourier transform of correlation 
functions should go like 



dN 



1 



The propagator in the above equation is 
defined as 



d 2 k t fcf + asu" 1 



(29) 



which plateaus off to a constant for small 
D(x + ,k,y + ,q) = f d 3 xd 3 y e~ ikx+iqy G(x, y) ,(24)nomentum. This is shown in Fig. 1. 

' Hence, asH is like a Debye scale which 

where guarantees the finiteness of the gluon dis- 

_ „ tribution functions for small momentum. 

G^( Xi y) = -i<<f> a (x)<f> P (y)> ■ (25) 

Now, in the correlation functions of the 
classical field (A t (x t )A t (y t )), the dejDfeW/d 2 ^. 
dence on x~ is only through a step func- 
tion, and upon Fourier transforming gives 
only a factor of l/k + . The distribution 
functions associated with this field, to all 
orders in a|/i 2 , have the general form [IT] 



dN 



7T i? 2 dxd 2 k t 



!:v;2 l)l --H(e t /a^)m 



7P 



x as 



When the relevant momentum scale is 
k 2 >> a|/i 2 , the theory is in the weak 
coupling region and may be evaluated per- 
turbatively. When y — > oo, 




H{y) = l/y. 



(27) 



and we obtain for the distribution function 
the Weizsacker- Williams result 



1 dN a S/ u 2 (iV 2 - 1) 1 



Ti R 2 dxd 2 q t 
scaled by /i 2 



7T^ 



xql 



2 > 



(28) 



For small values of k 2 « a|// , we are 
in the strong coupling phase of the the- 
ory. In this phase of the theory, we expect 
that there should be no long range order. 



Figure 1. Transverse momentum distribution 
for Weizsacker- Williams fields 



Note that the gluon distribution func- 
tion in the weak coupling Weizsacker- 
Williams regime is additive in A. If the 
behaviour in the strong coupling regime 
would resemble that in Eq. (|29|) , the gluon 
distribution function would then be pro- 
portional to A 2 / 3 . If we integrate over all 
k 2 to some scale Q 2 , the integrated distri- 
bution would probably have a dependence 
which lies between A 2 ^ 3 and A. One must 
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keep in mind however, the fact that higher 
orders will likely modify the A dependence 
of the classical distribution functions. 

4. Small fluctuations and higher or- 
ders in perturbation theory 

In previous sections, we discussed only 
the classical solutions of the Yang-Mills 
equations in the presence of the valence 
quark source defined in Eq. flT4]). In this 
section, we will outline a procedure to 
systematically compute quantum correc- 
tions to our background field to all or- 
ders. This procedure is exactly equivalent 
to the familiar Dyson-Schwinger expan- 
sion in many body physics ]20 . 



We begin by considering small fluctua- 
tions around our classical background, 

A cl = A° d + 6A. (30) 

Substituting this in the partition func- 
tion in Eq. (|I~2"D, we keep only terms upto 
0(5 A 2 ) in the action. The small fluctu- 
ations propagator may be computed di- 
rectly from the action or equivalently by 
substituting Eq. ([Kj) in the Yang-Mills 
equations, keeping terms linear in 5 A, and 
solving the resulting eigenvalue equation. 

The gluon Green's function is then ob- 
tained from the general expression 

SA x (x)SA* x (y) 



G(x,y) 



E 

A 



A 



(31) 



In Ref. (13), we obtained Green's func- 
tions for scalars, vectors and fermions 
in the non-Abelian Weizsacker- Williams 
(NAWW) background field. It turns out 
that the gluon Green's function we com- 
puted there in A + = gauge, is incor- 
rect and is actually the Green's function 



in A~ = gauge instead. The direct com- 
putation of the Green's function in Light 
Cone gauge A + = 0, is hampered by the 
presence of the singular source term in 
the small fluctuation equations of motion. 
The correct procedure for computing the 
Light Cone gauge Green's functions is dis- 
cussed in Ref. (14). The reader is referred 
to the papers Ref. (13) and Ref. (14) for 
relevant details. 

The general expression for the two point 
function to all orders in the coupling con- 
stant is given by 

(A t A t ) = (Af) (Af) + (A q t u A q t u ) , (32) 

where < Af > is the classical field to all 
orders and < Af ' A\ u > is the small fluctu- 
ation Green's function in the background 
field. The above expectation values also 
include a stochastic averaging over the 
sources of external charge with the Gaus- 
sian weight discussed in Eq. ([5|). Recall 
that as a consequence, there are two ex- 
pansion parameters in the theory: as and 
ag/i. So one should be wary of naive as 
counting! 

The zeroth order contribution to < 
Af > is the classical background field we 
discussed in previous sections and is of or- 
der 0(1/(7) J2i( a sf 1 ) n - The lowest order 
contribution to the distribution functions 
then gives 



(0) AO) 



At A 



oc 



g z 



■(a S fji) 



2/i 



asH 2 (n = 1) ,(33) 



which is the result obtained in Eq. (f28[). 

The next order contribution 0(g°) to 
the distribution function has two compo- 
nents. The first is from the small fluctua- 
tion Green's function (Af u Af u ) in Eq. (|32j) 
above. Recall that the Green's functions 
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we compute are directly related to distri- 
bution functions by Eq. (^). In momen- 
tum space this is equivalent to both legs of 
the distribution function having the same 
momentum. Performing the computation, 
we find a term proportional to log(x)/x. 
If we iterate this procedure to all orders 
(the distribution functions now forming a 
ladder)- summing the leading logs in x, 
we might get a series which looks like 

1 / C 2 a 2 \ 

- 1 + Ca s logx + — log 2 x + • • • = - 
x \ 2 J 

The crucial point is: what is the coefficient 
C? How does it depend on /z 2 ? This still 
remains to be settled: we are working on 
it! 

The other 0(g°) contribution to the dis- 
tribution function comes from the term 



as a function of fi 2 . This result is not en- 
tirely surprising 0] but it is nevertheless 
interesting to see how it arises in our for- 
malism. 

5. From virtual dream to hard real- 
ity 

In the preceding sections, we investi- 
gated the properties of the NAWW field 
1 of a^mgle, large nucleus. This is relevant 



A^AP in Eq. (0). Here Af includes 
the one loop radiative correction to the 
background field. One can show |T3[ from 
the operator equations of motion that 
satisfies the equation 

<9 2 A-' a + A(<M 1),a )+s/< 



abc 



lim d„ 



G%(x, 



The final term on the left hand side is the 
induced current in the background field. 
We find that it has exactly the same struc- 
ture as the classical current j£ and its 
only effect is to renormalize the exter- 
nal charge p(x t ). It therefore also follows 
that A^ has exactly the same structure 
as A[ \ The only change is that 

as — * ot s {^ 2 ) ■ 



(36) 



We may therefore conclude that the only 
effect of the classical field (Af) is to re- 
place the bare coupling constant in the lad- 
der sum by a coupling constant which runs 



x 1+c fer aeeply inelastic scattering off nuclei at 
small x. To understand nuclear collisions, 
however, we need to understand how the 
NAWW fields evolve after the collision. In 
addition, we need to address the question 
of how the virtual partons go "on shell" 
and undergo hard scattering. The scat- 
tering is hard because most of the wee 
partons are in the "hard" kinematic re- 
gion asfi C if < /i. In the first part of 
this section, we will discuss ideas of Larry 
McLerran and collaborators at Minnesota 
on the classical problem of the evolution of 
\ _thft-iNAWW fields. The latter part will be 
a brief discussion of "onium-onium" scat- 
tering as formulated by A. H. Mueller. 

5.1. The dynamical evo- 

lution of classical non— Abelian 
Weizsacker— Williams fields in 
nuclear collisions 

Recently, A. Kovner, L. McLerran and 
H. Weigert |22] managed to solve the clas- 
sical problem of the evolution of these 
fields. They showed that at late times, the 
equations for the evolution of the fields are 
linear. 

Below, we outline the important re- 
sults in their paper and refer the inter- 
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ested reader to their paper for further de- 
tails. The Yang-Mills equation for the two 
source problem is 



DpFi™ = J v {x) 



(37) 



where 

J± = 
J 1 = 



0. 



boost invariant variable r and has no de- 
pendence on the space-time rapidity vari- 
able y = -In 2=. That the solution is in- 

<? 2 x 

dependent of y suggests that the parton 
distributions will be boost invariant for all 
later times. This result therefore justifies 
Bjorken's ansatz p3 for the subsequent 



(38) 



In the above, pi(x t ) and P2(x t ) are the 
color charge distributions of the valence 
quarks of the two nuclei. 

Before the two nuclei collide (for times 
t < 0), the above equations of motion are 
satisfied by the classical field configuration 
where = and 



hydrodynamic evolution of the system. 

The equations that result for a{j, x±) 
and 03(1", x±) are highly non-linear. The 
detailed expressions are given in Ref.(22). 

Asymptotically, for large r, 



a(T,x±) 



V{x ± )e(r,x ± )V\x ± ) 



V(x ± ) 



e l (r,x ± )--d i 



where the value of the gauge transfor- 
A l = 6(x-)6(-x + )a\(x ± ) + 6(x + )9(-x~)a 2 (x L }^\ on is determined by the field 

equations. It results from solving the non- 
linear time evolution equations for the 
fields. 

The equations of motion for the fields 
in the asymptotic region are linear for e 
and e\ The solutions to these equations 
at asymptotically large r are of the form 



The two dimensional vector potentials are 
pure gauges (as in the single nucleus prob- 
lem!) and for t < solve 



v • 0:1,2 = g 2 pi,2(x±) 



(40) 



The interesting aspect of this solution is 
that the classical field configuration does 
not evolve in time for t < 0! This is a con- 
sequence of the highly coherent character 
of the wee parton cloud in the nuclei. 

The above solution for t < is a fairly 
straightforward deduction from the single 
nucleus case. What is very interesting is 
that the above mentioned authors find a 
solution to the field equations for t > 0. 
It is given by 



« a (r, x±) 



a a ' l (r, x±) 



d 2 k. 



1 



(27T) 2 
d 2 k ± 

WT 2 



2u 
1 



T 



3/2 



a a 2 {k A 



Jk ± -x ± - 



Ak 1 



T 



1/2 



In this equation, the frequency to 
and the vector 



K 



k j /u 



k 



(44) 



A ± 

A 1 

where r 



±x ± a(r, x_|_) 
al(r,x±) , 



V2 



(41) 



This 



solution only depends on the longitudinal 



The notation +C.C. means to add in the 
complex conjugate piece. 

With the above form for the fields, the 
expressions for the number and energy 
densities of the partons is straightforward. 
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For late times, near z = 0, the they look 



like 24 



dE 



dz 

T 



d 2 kj 
W) 



-UJ 



E I 



(45) 



8,6 



Recalling that ciy = <iz/£, we find that 
dE 1 



cfyd 2 fc ± (2vr) 3 



E I 



(46) 



Lb 



and the multiplicity distribution of gluons 
is 



dN 



dE 



dyd 2 kj_ uj dyd?k\ 



(47) 



The characteristic time scale for the dissi- 
pation of the non-linearities in the equa- 
tions for the time dependent Weizsacker- 
Williams fields can be estimated from di- 
mensional arguments. This is of the order 
r = 1/px = l/as/i. Note that for large 
/i's, i.e., large parton densities, this char- 
acteristic time gets smaller as one may in- 
tuitively expect. 

5.2. Onium— onium scattering 

An interesting re-formulation of the low 
x problem has been developed recently in 
a paper by A. H. Mueller ||25|| . He consid- 
ers an "Onium" (heavy quark-anti-quark) 
state of mass M for which as{M) -C 1. In 
weak coupling, the n- gluon component 
of the onium wavefunction obeys an inte- 
gral equation whose kernel in the leading 
logarithmic and large N c limit is precisely 
the BFKL kernel we discussed in the in- 
troduction. The derivation relies on a pic- 
ture in which the onium state produces 
a cascade of soft gluons strongly ordered 
in their longitudinal momentum; the i-th 



emitted gluon has a longitudinal momen- 
tum much smaller than the i — 1-th. 

The procedure to compute the wave 
function can be recast in terms of hamil- 
tonian perturbation theory in the follow- 
ing way: The Light Cone QCD Hamil- 
can be split into two pieces: 
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P n + P int , where P n and P f 



int 



toman 

P- = 

are the non-interacting and interacting 
pieces. Detailed expressions are given in 
Ref. (15). 

In the Infinite Momentum frame, where 
the parton picture can be applied, the in- 
teraction Hamiltonian acting on an onium 
state produces an state with a quark- 
antiquark pair and an extra gluon. First 
order perturbation theory can be used to 
give an expression for that state according 
to 



|0 



(i) 



E 

91.92 



< qi, <?2, p - qi - q2 P in t fo,P ~ h > 



(2 -part) 



P (0) ( a 
-^91.92 \y 



part) 



where E^(n — part) is the eigenvalue 
of P ~ on a state consisting of n free 
(quasireal) particles on mass shell; k\, 
p — h\ are the quark and antiquark mo- 
menta and the sum is over all gi, q2 con- 
sistent with overall momentum conserva- 
tion. Notice that the energy denominator 
in equation fl48|) will be equal to the (Light 
Cone) energy of the created gluon on mass 
shell, namely k 2 t /2kf, due to the strong 
ordering assumption. The state with n 
gluons and/or m quarks (anti-quarks) can 
be obtained from the vacuum by apply- 
ing n and/or m times the corresponding 
creation operators appearing in the Light 
Cone definition of the fields in the usual 
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way, 



d 3 k 



(2vr) 3 \/2fc 
d 3 k 



£e,(A,A;) \a a (X 



A=l 



(27T) 



b{k)e ikx + d\k)t 



-ikx 



The wave function in momentum space is 
obtained by projecting the state fl48|) onto 
the particular set of particle momenta. 
The gluon number, is obtained by squar- 
ing the wave function and adding its terms 
coherently. States with more than one 
gluon can be created by successively ap- 
plying the same procedure to the already 
created states, and the gluon density to 
that order can be computed by squaring 
coherently the wave function so obtained. 

In the large N c limit, the n gluons can 
be represented collection of n-dipoles. 
Hence, in high energy onium-onium scat- 
tering, the cross section is proportional to 
the product of the number of dipoles in 
each onium state times the dipole-dipole 
scattering cross section [p6| . This cross 
section is given by two gluon exchange || 
(the pomeron). More complicated ex- 
changes are also possible. 

It will be interesting to see if this ap- 
proach can be extended to the problem 
of nuclear scattering-how does it relate to 
the approach in the first part of this sec- 
tion? That there are parallels between the 
two approaches is evident-even though 
the large N c limit is a great simplification 
in the onium case. 

6. Summary and outlook 

In this work, we have discussed a many 
body approach to computing parton dis- 
tribution functions for large nuclei for x <C 



A^ 1 / 3 . Knowledge of these structure func- 
tions is crucial in formulating the initial 
fc)e*ten^itte(p&, fo)reultearelativistic nuclear col- 
lisions. Furthermore, such a first princi- 
ples calculation would eliminate the large 
uncertainities (especially at LHC) in jet 
cross sections and like rates f27fl . One 
may also apply these methods to under- 
stand the systematics of nuclear shadow- 
ing in deeply inelastic scattering experi- 
ments. Much work needs to be done in 
these directions. 

There are several open theoretical ques- 
tions which need to be addressed. A press- 
ing question is with regard to the applica- 
bility of weak coupling methods to large 
nuclei. Is fi 2 ~ A 1 / 3 large enough? For 
the largest nuclei, we have at most /i 2 ~ 2 
GeV 2 . At first blush, this seems extremely 
unlikely. However, due to the possibility 
of a hierarchy of scales in the problem (in 
analogy to finite temperature QCD), such 
a conclusion may be unduly pessimistic. 
If weak coupling methods are not applica- 
ble for computing the parton distributions 
of the largest nuclei, it is unlikely that 
they will ever apply during the subsequent 
stages of the evolution. This would im- 
ply that a weakly interacting Quark Gluon 
Plasma would never be formed! 

Another question we would like to ad- 
dress is whether the methods discussed for 
large nuclei can be extended to hadronic 
collisions. This depends on whether we 
can find a way to perform the quantum 
mechanical sum over the external sources 
in the partition function in Eq. (|8|). 

The x — > limit of QCD is a very inter- 
esting one because many features of the 
theory simplify in this limit. For instance, 
it has been shown recently that this limit, 
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for N, 
ory 



oo is an exactly solvable the- 



2S 



It will interesting to see whether 
these results can be recovered in the many 
body formalism discussed in this paper. 
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